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1 Introduction 

Our program of research aims at developing a stochas- 
tic description of the residual acceleration field onboard 
spacecraft (5-jitter) [1, 2] to describe in quantitative de- 
tail its effect on fluid motion [3, 4, 5]. Our main premise 
is that such a statistical description is necessary in those 
cases in which the characteristic time scales of the pro- 
cess under investigation are long compared with the cor- 
relation time of ^-jitter. Although a clear separation 
between time scales makes this approach feasible, there 
remain several difficulties of practical nature: (i), g- 
jitter time series arc not statistically stationary but rather 
show definite dependences on factors such as active or 
rest crew periods; (ii), it is very difficult to extract re- 
liably the low frequency range of the power spectrum 
of the acceleration field. This range controls the magni- 
tude of diffusive processes; and (iii), models used to date 
are Gaussian, but there is evidence that large amplitude 
disturbances occur much more frequently than a Gaus- 
sian distribution would predict. The lack of stationarity 
does not constitute a severe limitation in practice, since 
the intensity of the stochastic components changes very 
slowly during space missions (perhaps over times of the 
order of hours). A separate analysis of large amplitude 
disturbances has not been undertaken yet, but it does not 
seem difficult a priori to devise models that may describe 
this range better than a Gaussian distribution. The ef- 
fect of low frequency components, on the other hand, is 
more difficult to ascertain, partly due to the difficulty as- 
sociated with measuring them, and partly because they 
may be indistinguishable from slowly changing aver- 
ages. This latter effect is further complicated by the 
lack of statistical stationarity of the time series. 

Recent work has focused on the effect of stochastic 
modulation on the onset of oscillatory instabilities [6] 
as an example of resonant interaction between the driv- 
ing acceleration and normal modes of the system, and 
on cavity flow [7] as an example of how an oscillatory 
response under periodic driving becomes diffusive if the 
forcing is random instead. This paper describes three 
different topics that illustrate behavior that is peculiar to 


a stochastic acceleration field. In the first case, we show 
that g-jitter can induce effective attractive or repulsive 
forces between a pair of spherical particles that arc sus- 
pended in an incompressible fluid of different density 
provided that the momentum diffusion length is larger 
than the inter particle separation (as in the case in most 
colloidal suspensions). Second, a stochastic modulation 
of the control parameter in the vicinity of a pitchfork or 
supercritical bifurcation is known not to affect the loca- 
tion of the threshold. We show, however, that resonance 
between the modulation and linearly stable modes close 
to onset can lead to a shift in threshold. Finally, we 
discuss the classical problem of vorticity diffusion away 
from a plane boundary that is being vibrated along its 
own plane. Periodic motion with zero average vorticity 
production results in an exponential decay of the vortic- 
ity away from the boundary. Random vibration, on the 
other hand, results in power law decay away from the 
boundary even if vorticity production averages to zero. 

2 Acceleration induced interactions between 
pairs of particles 

Consider an ensemble of spherical particles of radii Ri 
and density p p suspended in an incompressible fluid of 
density pf and shear viscosity p. If the fluid is enclosed 
by perfectly rigid boundaries, the buoyancy force acting 
on each suspended particle is F = | n(p p — pf) R 3 g(t), 
where g(t) is the effective acceleration field. In the frame 
of reference co-moving with the container enclosing the 
fluid, g(t) is a body force, with intensity equal to the 
value of the acceleration of the container. For containers 
of reasonable size in a microgravity environment, g can 
be assumed to be spatially uniform. The hydrodynamic 
interaction between two such particles is given in the 
overdamped limit of Stokes flow by [8, 9], 

(It — * — * 
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where f is the relative position of particle 2 with respect 
to particle 1, Fi is the force acting on the i-th particle, 
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3 STOCHASTIC RESONANCE AND BIFURCATIONS 


and u >ij arc hydrodynamic mobility tensors, given, e.g., 
in references [8, 9]. After some straightforward algebra, 
the leading contribution at distances large compared to 
the particle radii is given by, 


(Pp ~ Pj) («? ~ fij) 



•m> (2) 


where 1 is the identity tensor. The first term in the right 
hand side describes the relative motion of two indepen- 
dent particles of different size. Both the longitudinal and 
transverse components of the second term in the right 
hand side of Eq. (2) are of the form, 


where, for the longitudinal component, A = (p P -pf)(Ri- 

R\)l 3 m * 

Consider an initial interparticlc separation r 0 > Ri. 
In this case, and for times shorter than the average time 
needed for the two particles to coalesce, the quantity y = 
r 2 /2A is a Wiener process if g(t) is Gaussian and white 
(( g } = 0, { g{t)g(t') ) = 2 DS(t - t r )\ and therefore the 
conditional probability for r is, 


P(r,t\ro,to) 


2\ 2 


(— o 2 ) 


\A\y/4nD{t - to) 


16DA 2 (t-t 0 ) 


(4) 


and the fact that the noise is Gaussian and white, one 
finds, 


dt - AD \ 8g(t) / ’ 


where S/Sg[t) stands for functional derivative with re- 
spect to g. Directly from Eq. (3), wefindthat<5(l/r(t))/<Sp(t) 
-A/r 3 , and therefore. 


d{r) 

dt 



( 8 ) 


identical to Eq. (6) with 1/t 3 replaced by { 1 jr 3 ). Since 
r is a positive quantity, d ( r ) /dt < 0 for all values of r. 
It is also interesting to note that the effective attractive 
interaction is not confined to the term proportional to 
1 jr in the hydrodynamic mobility, but that attractive 
contributions arise from higher powers of 1/r as well. 
In fact, this attraction is generic for over damped motion 
and multiplicative noise provided that the mobility is a 
decaying function of the interparticle separation [12]. 

The question naturally arises as to the behavior of 
pairs of particles near contact, or of particles near a 
solid wall. In either case, lubrication theory allows 
the calculation of the mobility tensor. The longitudi- 
nal component vanishes linearly with interparticle dis- 
tance whereas the transverse component becomes non- 
analytic (diverges logarithmically at short distances) [9]. 
In both cases, the mobility increases with interparticle 
separation leading to an average repulsion (d(r) /dt > 
0) following the same arguments given above. 


The ensemble average of r, ( r ) can be computed ana- 
lytically, 3 Stochastic resonance and bifurcations 


(r ) = 


■ [2D(t - to)] e 


4 

Jji 
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Consider the normal form appropriate for a pitchfork 
bifurcation in which A is the linearly unstable mode and 
B is some linearly stable mode. Then, 
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where D p (z) is a parabolic cylinder function [10] (for- 
mula 9.240). For short times, the asymptotic form of 
D p (z) for large 2 allows the computation of ( r ) 


mu 

77712 

r a ~ 

77121 

77722 

[ B 


9{t), 


+ 

(9) 


, , /\ A 2 D(t-t 0 )\ 

( r ) = 7*0 I 1 -4 I, ( 6 ) 

which decreases in time regardless of the sign of A. 
Therefore g-jitter induces an effective hydrodynamic at- 
traction between pairs of particles. The attractive inter- 
action is not confined to short times, but it arises directly 
from the 1 jr dependence in Eq. (3). By taking the aver- 
age of Eq. (3), using the Furutsu-Novikov theorem [11] 


with c 1 and the remaining coefficients of order one. 
We further assume that D ^ 0(e). Close to threshold, 
A changes in a slow scale T = et, compared with either 
the relaxation of B (of order one) or the noise. 

If the coupling to the unstable mode B is neglected, 
it is known that the threshold in the stochastic case re- 
mains at e = 0 [13, 14]. This can be seen by averaging 
Eq. (9) over the fast time scale so that (A(T)g(t)) ^ 
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4 FLOW DUE TO A RANDOMLY VIBRATING PLANE BOUNDARY 


A(T)(g{t)) = 0. The resulting equation for A(T) no 
longer depends on the noise. If, on the other hand, both 
equations are averaged over the fast time scale, we find 
that, 


(Bg) = (B)(g) + D{^J = Dm 2 iA + ... (10) 

where we have used the Furutsu-Novikov theorem [15, 
16]. Therefore, the correlation of B(t)g(t) itself evolves 
over the slow time scale T as a consequence of the 
fact that the equations for both A and B contain ex- 
actly the same stochastic process. The coefficient of 
the linear term in the equation for A(T) is now e + 
Dm 2 im \2 and hence the bifurcation point will occur at 
€ = -Dm.\ 2 TTl 2 l < 0. 

4 Flow due to a randomly vibrating plane 
boundary 

Consider a semi-infinite fluid layer occupying the re- 
gion x > 0 and a solid boundary at x = 0 which is 
being displaced along the y direction with a prescribed, 
time-dependent velocity vo(t). For an incompressible, 
Newtonian fluid, the y component of the velocity field 
in the fluid v satisfies, 


dv _ d 2 v 
dt V dx 2 ' 


( 11 ) 


an approximate description of the spectral components 
of the residual acceleration field in microgravity, as well 
as a convenient way to interpolate between the white 
noise and monochromatic noise limits. Narrow band 
noise is a Gaussian process that satisfies, [17] 

M*)> = o, (t>o(t)«o(0) = t>o e ~ l<- *' l/T cosfl(< - t'), 

(14) 

where ( ) denotes an ensemble average. Q is a character- 
istic angular frequency that corresponds to the peak in 
the spectral density of the process, and r is a correlation 
or coherence time determining the width of that peak. In 
essence, this process describes a periodic signal of char- 
acteristic frequency Q, but that only remains coherent 
for a time r on average. The amplitudes are distributed 
gaussianly with variance vq. The white noise limit is 
obtained when Dr — > 0 while vqt — D remains finite, 
whereas the monochromatic noise limit corresponds to 
Hr -* oo with Vq finite. We do not attempt to find a 
general solution of Eq. (1 1) for an initially quiescent 
fluid and narrow band forcing. We focus instead on long 
time or statistically stationary averages. We find, 

/ °° faM 

du}P(uj)e v * . (15) 

- oc 

where P{uj) is the power spectrum that corresponds to 
the autocorrelation function (14). The integral can be 
carried out explicitly in the limits of small and large r. 
In the vicinity of the white noise limit, we find, 


with boundary conditions v(x = 0, t) = uo(f) and i>(z -* 
oo, t) = 0. We consider a function u 0 (t) which is a 
stochastic process in time, and hence proceed to solve 
the initial value problem (11) with a stochastic boundary 
condition. Let v Q (t) be a Gaussian, white stochastic pro- 
cess with mean (v 0 (t)} = 0 and variance (vo(f)uo(O) = 
2 DS(t - t'). We find that, 


(r 2 (x,<)) 


2 Dv 
nx 2 



( 12 ) 


{v 2 (x)) = 


2isvlr 1 i>o(120 -I- 4Q 2 x 4 /f/ 2 )^ 3 


- 0 ( 0 . 
(16) 

The first term in the right hand side is the white noise 
limit already given in Eq. ( 1 3). The first correction term 
is also a power law decaying as x -6 away from the wall. 
The low frequency part of the power spectrum dominates 
the decay of the velocity field at long distances and leads 
to a very slow rate of decay. In the opposite limit of 
rrightarrowoc, we find, 


The velocity disturbance propagates into the fluid dif- 
fusively, with a diffusion coefficient v. At long times, 
however, even though the vorticity produced at the wall 
averages to zero, the velocity does not decay exponen- 
tially away from the wall, but rather as a power law, 

(v 2 (x,t -+ oo)) = (13) 

We now turn to the case in which the velocity of the 
boundary is not white, but a narrow band noise instead. 
As discussed elsewhere [2], narrow band noise provides 


(v 2 (x)) = («o) e + 



e _ Iy /m±uzi 

\ 



e V ^ | 


In the limit r -> oo the term within braces vanishes and 
one recovers the classical result that the characteristic 
velocity \J(v 2 ) decays exponentially into the fluid with 
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a characteristic decay length given by If, on the 

other hand, the spectrum of the forcing function has a 
finite width, the velocity field still d ecays exponen tially, 
but the decay length increases to yj2vjjQ. - 1/r). 
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